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We find analytic continuous wave (cw) solutions for spinor Bose-Einstein condensates (BECs) in 
... a magnetic field that are more general than those published to date. For particles with spin F = 1 

CO ' in a homogeneous one-dimensional trap, there exist cw states in which the chemical potential and 

, wavevectors of the different spin components are different from each other. We include linear and 

■ quadratic Zeeman splitting. Linear Zeeman splitting, if the magnetic field is constant and uniform. 
(N . can be mathematically eliminated by a gauge transformation, but quadratic Zeeman effects modify 

I the cw solutions in a way similar to non-zero differences in the wavenumbers between the different 

. spin states. The solutions are stable fixed points within the continuous wave framework, and the 

' coherent spin mixing frequencies are obtained. 

^ ' PACS numbers: 03.75.Hh, 03.75.Kk, 03.75.Mn 

^ ! I. INTRODUCTION 

. . . 

^ I ^ I Atomic Bose-Einstein condensates (BECs) with nonvanishing total angular momentum have been experimentally 
' obtained, in which the BECs contain several spin components, i.e., they are spinor BECs. Examples are H (with 
d , total spin F = and F = 1) ''Li (F = 1, F = 2; note that the scattering length is negative, making the BEC 
^ • self-attractive) d ^^Na {F = 1, F = 2) H], ^^K (F = 1, F = 2; scattering length has a small negative value) [H, 
^ 2, F = 2) 0, s^Cr (F = 3; with a large magnetic moment) 0, i], ^^Rb (F = 2, F = 3; negative scattering 
■ ~ length), 87Rb (F = 1, F = 2) ^^^Cs (F = 3, F = 4) [H, ^^^By (F = 8) [l|, and ^^^^Er (F = 6) [HI. When a 

BEC is held in a magnetic trap, only spin components of one sign arc trapped, often leaving just one spin orientation. 
^ • An optical trap, however, may hold all spin components of a given hyperfine state. In this case, the spinor character 
1^ I is important, and a scalar model of a BEC is insufficient. Optical traps for spinor BECs are now common, and there 
. is work on creating BECs using only optical traps 0, Q . Optical traps have also been created in the form of a closed 
O ; ring [IHp. 

, ^ , , The spinor properties of BECs have been the focus of recent theoretical and experimental research [2l|, [g^] , especially 
^^Na, *^Rb, and ^^Cr. The spinor properties of a BEC can critically affect the dynamics: The spinor character of 
CO i BECs underlies the formation of domain walls (transition regions, which may be stable or unstable, between distinct 
K*" ' spin domains) and vortices [ll|, HI, [l^j- Oscillatory coherent spin mixing occurs in spinor BECs [25|, [2^. Spinor 
BECs are subject to modulational (Benjamin-Feir) instabilities [2^ |27h29| | even when the nonlinearity is repulsive, 
whereas scalar BECs are not. Spinor BECs support a variety of soliton solutions that are not found in scalar BECs 

[MH,!!!!!!. 

Here we study continuous wave (cw) [plane wave] solutions of F = 1 spinor BEC condensates. Section HIl introduces 
the model, which includes spin-dependent and spin-independent mean-field effects, and linear and quadratic Zeeman 
I effects (but docs not include spin-dipolar effects, which can be important, e.g., in ^^Cr). Section Hill derives the most 
general possible cw solutions for F = 1 spinor BECs on a homogeneous background with a homogeneous magnetic 
field. This section gives, even before inclusion of the linear and quadratic Zeeman effects, more general families of cw 
solutions than have heretofore been published. Section [IVI contains a summary and conclusions. 

>< 

■ II. QUANTITATIVE MODEL FOR SPINOR BECS WITH MAGNETIC FIELDS 
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The Hamiltonian density for an F = 1 spinor BEC with linear and quadratic Zeeman effect (and without significant 
spin-dipolar coupling) is [sgI . |37| 

n = |^V<i>t . V<i>a + ^|$t$t$,ci,^ + ^^l^lFat ■ F,,5,$,,$, _ pB<^>iFl,<^>, + qB''^i{Fl,)H,, (1) 

where $ = (0i,0Oi0-i)' is a vector with the amplitudes of the Mp = 1 component (spin in the same direction 
as the magnetic field). Alp = 0, and Mp ~ —1; m is the mass of the atom, cq and C2 are the coefficients of the 
spin-independent and spin-dependent parts of the mean-field; cq gives rise to self-phase modulation and C2 is the 
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spin-dependent mean-field coefficient wfiicli also gives rise to phase modulation, and parametric nonlinearity. 
vector in which each component is a 3x3 spin-1 matrix, i.e., the dimensionless spin vector F has components. 






= -= 1 1 , Fj, = -= i , = . (2) 



The magnetic field B is taken to be constant and uniform and in the z-direction, B = Bz, and p and q are linear 
and quadratic Zeeman coefficients [23l - f25| . If the BEC is confined to one dimension (by a strong optical trap in the 
transverse directions), the Hamiltonian gives the following governing equations for the BEC: 

"-^wA = -T^TT^^i+cod'/'ip+i'/'oi'+i'/'-in^i 

at 2m oz^ ^ ' 

+C2 [(|</'i|' + |0oP - l^-iP) 4>i + '^oCi] + {-VB + 952)01 , (3a) 



+C2 + |0-iP) 00 + 20i0;^0_i] , (3b) 



+C2 [(-I'/'il' + |0ol' + (^-1 + <l>lri] + {pB + qB^)cjy_^ . (3c) 

Time and space are denoted by t and z, respectively. The system has Galilean invariance. Materials with negative C2 
are ferromagnetic, since with C2 < 0, at a given particle density and at zero magnetic field, the energy density is lower 
when the local BEC is in a pure spin state Mp = 1 or Mp = — 1, than in a state composed of a mixture of spins Mp = 1 
and Mp = — 1. The BEC prefers to have a non-zero net spin locally. Materials with positive C2 are antiferromagnetic, 
or polar, since with C2 < 0, at a given particle density and without a magnetic field, the energy density is lower when 
the local BEC is 50% Mp = 1 and 50% Mp = — 1 . The BEC prefers to be in a state where opposite spins balance each 
other exactly and, cancel out in the total magnetic moment. The nonlinear coefficients, which are proportional to the 
s-wave scattering lengths oq, 0,2 for the total spin F = and F ~ 2 channels. 170 — {'^'kTi^ /m)aQ, 32 = (47r?i2/m)a2, 
with the nonlinear coefficients in the governing equations above given by cq = (go + 2(72)/3, C2 = —{go — 52)/3. If the 
BEC is confined to one dimension, the values of the nonlinear coefficients are modified as discussed in Refs. [s^. [39l | . 
For *''Rb, the scatterings lengths are oq = 101.8 as and 02 = ao — 1.45 as, where as is the Bohr radius [40l442| . This 
gives ®''Rb a negative C2, which makes it ferromagnetic. For ^'^Na, the scattering lengths have been measured to be 
ao = 50.0 as, a2 = ao + 5.0 as. ^"^Na has a positive C2, so it is antiferromagnetic, or "polar." The masses of the 
two atoms are mRb 1.4192 x lO'^^g, mNa = 3.817 x lO'^^g. In bulk, these yield co(Rb) = 5.26 x 10~^^ergcm^, 
C2(Rb) = -2.52 X 10~^°ergcm3 and co(Na) = 10.33 x lO'^*^ ergcm^, c2(Na) = 3.23 x lO'^^crgcm^; the rations C2/C0 
are -0.0048 for *^^Rb and 0.0313 for ^SNa. The quadratic Zeeman coefficient is g = ft. x 575Hz/G2 for *^Rb, and 
g = /i X TOHz/G^ for ^'^Na. Equations ^ are integrable when C2 = (in which case the system is a set of generalized 
Manakov equations [H, Hil ) or C2 = co PO, 31, 45||. 

The linear Zeeman splitting can be eliminated from the governing equations ([3]) by the change of variables (i.e., 
the gauge transformation) <f>i = ipi cxp[i(pB /h)t], <j)Q = i/'oi 4'-i = V"-! exp[— z(p_B/?i)t]. Let us take this as done, but 
retain the same variable names as previously, to avoid excessive complicated notation. Because the linear Zeeman 
splitting can be eliminated by this gauge transformation, the analysis of the nil linear Zeeman splitting, p = 0, can 
equally well describe the case with non-zero linear Zeeman splitting. Generally, analyses should either make use of 
this gauge transformation or not force all the spin components to have the same frequency (chemical potential), or 
one may erroneously impose restrictions that do not exist in the physics. Quadratic Zeeman splitting cannot be 
eliminated by a change of variables. 

The analysis below is in terms of infinite continuous (plane) waves. It applies equally to waves with periodic 
boundary conditions, either physical ones due to the BECs existing in a circular trap [l6|-|20|, or theoretical ones due 
to the need to limit the numerical domain. Periodic boundary conditions quantize the wavenumbers. 

The governing equations ^ can be non-dimcnsionalized by the change of variables 

t' = t/U, (4a) 
z' = z/zd = z/y/htd/m, (4b) 
= 0j70d = 4>iN'hl[cQtd). (4c) 

In dimensionless variables (|lcl) . the governing equations (jS]) take the same form but with ?i— !-l,m— )-l, co— s-1, and 
C2 — ^ C2/ cq. The physical frequencies and wavenumbers are equal to the dimensionless quantities times t^^ and = 
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{yjhtd/m) ^, respectively. The physical amplitudes of the BEC are the dimensionless ones times yjfi/ (cotd)- The 
physical energy is equal to the dimensionless energy multiplied by h/td, and the physical energy is the dimensionless 

quantity times hT^^'^nv'I'^iJ''^'^ . Wc will use the dimensionless variables in the figures in order to emphasize the 
generality, but retain the dimensions in the body of the text. 
The Lagrangian density that gives Eqs. ^ is 



- ^ + 0500,* - 0O0S,t + - <^-\<t)*-x.t) 

2\2 



+ |0or + |0-i|^ 



C2 



(5) 



and the Hamiltonian density is 



2m 2 



C2 



10- 



^'t0O<?^ll + 010S^0-1 



l^-iP). 



(6) 



If any two of the three spin fields are zero, then the remaining field is governed by a simple nonlinear Schrodinger 
equation, which is completely integrable [l^S^]- If the Mp = field is nil ((/)o = 0), then the Mp = ±1 fields (0i, 
0_i) are governed by a pair of coupled nonlinear Schrodinger equations, 



^Tf^'l'^ + [(^0 + C2)|01p + (co - C2)|0-lP] 01 , 
-^TTl^-^ + [(^0 - C2)|01p + (co + C2)|0_lP] 0_i 



(7a) 
(7b) 



The coupled nonlinear Schrodinger equations have been intensely studied (see, e.g., Ref. [43|). This case is completely 
integrable if and only if C2 = 0, cq — C2, or cq = (ioj . Depending on the coefficients, there may be bright soliton 
solutions [i^ [sO] , dark soliton solutions , bright-dark soliton solutions [H, HI] , and domain wall solutions [13, [11] . 
We will concentrate less on this limit, and more on aspects of spinor BECs that do not overlap with the thoroughly 
studied coupled nonlinear Schrodinger equations. 



III. CONTINUOUS WAVE SOLUTIONS 

Continuous waves (cws) are the simplest shapes, so cw solutions are where analysis of the system should begin. 
The most general cw is 

01 = Aicxp[i{9i + kiz ~ ujit)] , (8a) 

00 = Ao exp[i(0o + kaz - ujot)] , (8b) 

0-1 = A_i exp[i((?_i + fc_iz — w-i^)] , (8c) 

where the parameters are real- valued and, without loss of generality, Ai, Aq, A-i are positive definite. Note that this 



is more general than the cw ansatz in Ref. |3l| . in that the frequencies and wavenumbers need not all be the same. 



The analysis herein shows that there exist a wider range of cw solutions than, for example, in Ref. [31[ . Put the trial 
function © into Eqs. 

If C2 = 0, then Eqs. ^ are a generalized Manakov system [i^, but with three fields rather than two. This limiting 
case is completely integrable [4J]. There are cw solutions for every value of the amplitude (Aj), wavenumber (kj), 
and phase {9j ) . The frequencies of the fields are 



tip 

hu;, = ^+coiAl+Al+Al,). (9) 
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If C2 7^ 0, then the parametric term requires a relation between the phases of the three fields, 

1 



ko 



UJQ 



(fci +fc_i), 
(wi + W-l) , 



where n is an integer. The equations for the magnitudes of the fields are 



Tiuji = 
huJo = 



H 

2m 

H 

2m 

f^kl^ 

2m 



- CO (Aj + AI + Al,) + C2 [Al + Al- A\ + (-1)"A2A_iMi] + qB'A^ , 

- CO +Al^ A\) + C2 [A? + A\ + 2(-l)"^a_i] , 

+ CO {A\ + ^2 + ^2_^) ^ [_^2 ^ ^2 ^ ^2_^ ^ (-1)"A2AiM_i] + gB^A.i 



Equations PT|) and (jlObp give a formula for the magnitude Aq of mode (/>o, 



(?iV2m)[(fci ~fc_i)/2]2 + gB2 



= 2(-l)"AiA_i 1 



C2[Al + (-l)"A_i]2 



(10a) 
(10b) 
(10c) 

(11a) 
(lib) 
(11c) 

(12) 



For a cw solution to exist, Aq must be non-negative, as per ansatz ^|o);, bu itn nquaic munt uc uun-ncgauivc, ^o 
The ranges in which the different cw solutions exist depends on the magnitude and sign of the quantity \^(k 



, so its square must be non-negative, An > 0. 



k-rf 



qB ]/c2. If the quadratic Zeeman coefficient is non- negative {q > 0), the sign of this is either zero or 



equal to the sign of C2. If the quadratic Zeeman coefficient is negative (q < 0), a sufficiently strong magnetic field can 
change the sign of the quantity. For particles with F = I, the quadratic Zeeman coefficient is normally positive, q > 0. 
(Particles with F = 2, normally have negative quadratic Zeeman coefficients, q < 0.) The sign and magnitude of q can 
be changed, however, by using the alternating current Stark shift with microwave radiation [l2, 21, 26]. We will refer 



to [h (fci — /s_i)2/8m -|- qB'^]/c2 > as the generalized antiferromagnetic case, and [Ti (fci — k-if 



gB2]/c2 < 



as the generalized ferromagnetic case. In the generalized antiferromagnetic case, there are cw solutions with even n 
(which we will denote by n = in the labels below) when 



(^1 



A.,f>- 



1 

C2 



2m 



fci — fc_ 



+ qB' 



(13) 



and there are cw solutions with odd n (which we will denote by ri = 1) when 

1 



0< (Ai-A_i)2 < 



C2 



/ fci — fc_ ^ ^ 

2m 



qB^ 



(14) 



There are no n = 1 type cw solutions when Ai = A_i because the difference of the amplitudes in the denominator of 
Eq. ()12p makes the value of Aq go to infinity as [Ai — A-i) approaches zero. In the generalized ferromagnetic case, 
there are n = cw solutions for all values of the spin Mp = and spin Mp = —1 magnitudes {Ai, A_i £ TZ-^). The 
generalized ferromagnetic case does not support any n — \ type cw solutions. Recall that there exist cw solutions with 
vanishing Mp — BEC fields for any values of the Mp = ±1 BEC magnitudes, whether the BEC is ferromagnetic or 
antiferromagnetic. This existence ranges are summarized in Table H] and in Fig. [1] 
The Hamiltonian density of the cw solutions is 



n = ^{k\A\ + klAl^k\A\)^''-l[A\ 



A^_,f 



C2[i(A 



At 



Al{A\ + A\ + (-l)'MiA_i)] + qB^^Al + A\) , 



(15) 



where Aq may either be zero or one of the non-zero cw solutions [Eq. (|12p]. and fcp = (fci + fc_i)/2. 

Note that all three families of cw solutions include solutions in which the spin components have mutually different 
wavenumbers - the solution with identical wavenumbers is a limiting case. This implies that the different components 
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X = [h^{ki ~ k-if/8m + qB^]/c2 


cw (CNLS) 


cw (n = 0) 


cw (n = 1) 


> 




[Ai+A-if > X 


< (Ai - A-if < X 


< 


Ai,A-i G TZ^° 


Ar,A-i G 7^^'' 


No solutions. 



TABLE I: Existence ranges of tiie three possible cw solutions, for different values of the spin-dependent nonlinear coefficient 
C2, particle mass m, quadratic Zeeman coefficient q, magnetic field B, and given difference in the wavenumbers between the 
Mf = ±f spin components. The allowed magnitudes of the amplitudes of the fields for the BEC with spin Mf = ±f are 
denoted by A±i, as in the ansatz ((8]), and the magnitude of the amplitudes of the BEC with spin Mf = is either ("CNLS") 
or goes according to the formula (|f 21) . 




FIG. 1: (Color online) Existence ranges of the cw solutions in spinor BECs in the generalized antiferromagnetic case (Arcf = 
)^/8m + qB'^\/c2 > 0; note that in the most familiar cw solutions, with fci = fco = and zero quadratic Zeeman 
splitting, Arof = 0; in this limit, the n = cw solutions have Ao = ^2A\A-\ and exist for all values of the particle densities 
Mf = ±1, and the n — 1-type cw does not exist). Part (a) shows cws with relative phase of the Mf ~ field corresponding 
to even n in the trial function Eq. (|8]), with 8i + 9-i — 26o = mr. Part (b) shows the cws with relative phase designated by 
odd n. For odd n, the particle density of the AIf ~ field approaches infinity as the densities of A4f = 1 and A4f = — 1 fields 
approach equality. 

may move relative to one another. In the mathematical model, the domain is infinite. In an experiment with cws with 
different wavenumbers in the different components, one would have to either ensure that the fields remain in a specific 
location for long enough, either by doing the experiment before the edges encroach on the middle, replenishing the 
fields at the edges, or arranging the fields in a ring (confined by a toroidal potential) [l6l - [20| . 

Figures [2] through |6] show, for representative cw solutions, the values of the amplitudes of the spin Mp = fields 
and the Hamiltonian densities. In these sets of figures, the nonlinear coefficients are those of ^^Na and ^^Rb, and 
the the plots are scaled such that they represent any generalized antiferromagnetic and ferromagnetic cw solutions 
(though excluding the sign changes that can come from a negative quadratic Zeeman coefficients with a non-zero 
magnetic field). 

We carried out a stability analysis of the cw solutions within the continuous wave assumption, i.e., without looking 
for modulational (Benjamin- Feir) instabilities or sound waves (alias Bogoliubov excitations or phonons), that is, 
small perturbations of arbitrary wavenumber on top of the cw solutions |56H59f . A study of the sound waves and 
modulational stabilities is a large and lengthy topic. We will defer these results to another article. References (22l.[27l- 
[29j analyze sound waves and modulational instabilities on top of subsets of the complete set of possible cw solutions 
herein. The stability analysis generalizes the ansatz ([5]) by allowing the variables for the magnitudes of the amplitudes 
Ai^o,-ij their phases 0i.o,-i, and wavenumbers fci.o,-i to vary with time. Inserting this ansatz into the dynamical 
equations ([3]) gives the results that the total particle density, N = Af + Aq + A^j^, and the magnetization, M = 
A\ — A^_-^, are conserved. Moreover, the wavenumbers ki^^-i are do not change with time, due to the physics with 
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FIG. 2: (Color online) Hamiltonian density of cw solutions of a BEC of ^'^Na atoms in their ground state, having no particles 
with spin component AIf — 0. The quantities are shown in dimensionless units [see Eq. (|4cp ]. for generality and to avoid large 
amounts of notation. The cws in the figure have a fixed unit difference between the AIf = 1 and Mf = — 1 spin components 
wavenumbers, ki — = 1, and the magnetic field is zero. The Hamiltonian density is shown for a range of magnitudes of the 
BEC fields with Alp = ±1. The (dimensionless) reference amplitude on the axes is 




FIG. 3: (Color online) Amplitude of the spin Mf = component and Hamiltonian density for cw solutions of a BEC of 
^^Na BEC, with the relative phase in the spin AIf ~ component corresponding to even n in the trial function ((H]), with 
9i + 9-1 — 29o — mr. The quantities are shown in dimensionless units [see Eqs. H4c|) ]. The cws in the figure have a fixed unit 
difference between the Mf ~ 1 and Mf = — 1 spin components wavenumbers, fci — fc_i — 1, and the magnetic field is zero. 
The magnitude of the spin Mf = amplitude and the Hamiltonian density of the cw is shown for a range of magnitudes of 
the BEC fields with Mf = ±1. The (dimensionless) reference amplitude is Ard = [h^{ki — fc_i)^/8m + qB'^]/c2. 



the cw assumption. We may reduce the osciUations about the fixed points that represent the cw solutions to the 
dynamics of two variables, ^o(i)^ and [9i{t) + 0-i{t) — 29o{t)] which are conjugate to each other. The fixed points 



are stable to small perturbations, with frequencies 



/ ,2 



3(-l)" + 2 



A. 

^-1 



(16) 



The frequencies are always real- valued. Thus, all the n = 0, 1 cw solutions are stable within the continuous wave frame- 
work (i.e., not considering modulational instabilities, which due to length is deferred to another paper). Physically, 
these are the frequencies of oscillatory coherent spin mixing 25, 26, 33]. 

This oscillation is not the same as that in Refs. [U and[35|. In those papers, the frame of reference can be rotated 
such that the solutions have no BEC component with spin Mp = 0. The governing equations then reduce to two 
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FIG. 4: (Color online) Amplitude of the spin Mf = component and Hamiltonian density for cw solutions of a BEC of 
■^^Na BEC, with the relative phase in the spin Mf = component corresponding to odd n in the trial function Eq. (|8]), with 
Gi + 6-1 — 200 = nvr. The quantities are shown in dimensionless units [see Eq. (|4c|) ]. The cws in the figure have a fixed unit 
difference between the Mf = 1 and Mf = — 1 spin components wavenumbers, fci — fc-i = 1, and the magnetic field is zero. 
The magnitude of the spin Mf = amplitude and the Hamiltonian density of the cw is shown for a range of magnitudes of 
the BEC fields with Mf = ±1. The (dimensionless) reference amplitude is A^^i = — fc_i)^/8m + g_B^]/c2. For odd n, 

the particle density of the A'/f ~ BEC (as well as its Hamiltonian density) approaches infinity as the values of the densities 
of Mf = 1 and A'If = — 1 fields approach equality. 








FIG. 5: (Color online) Hamiltonian density of cw solutions of a BEC of *^Rb atoms in their ground state, having no particles 
with spin component Alp ~ 0. The quantities are shown in dimensionless units [see Eq. (|4c|) ]. for generality and to avoid large 
amounts of notation. The cws have a fixed unit difference between the Mf = 1 and Mf = —1 spin components wavenumbers, 
ki — k-i = 1, and the magnetic field is zero. The Hamiltonian density is shown for a range of magnitudes of the BEC fields 
with A/f = ±1. The (dimensionless) reference amplitude on the axes is Ard = —[^^{ki — k-i)^ /8m + qB^]/c2. 



8 




FIG. 6: (Color online) Amplitude of the spin Mf = component and Hamiltonian density for cw solutions of a BEC of 
*^Rb BEC, with the relative phase in the spin Mp = component corresponding to even n in the trial function Eq. ((8]), with 
9i + 6-1 — 29o = nn. The quantities are shown in dimensionless units [see Eq. H4c[) ]. These cws have a fixed unit difference 
between the Mf = 1 and AIf = —1 spin components wavenumbers, ki —k-i = 1, and the magnetic field is zero. The magnitude 
of the spin Mf = amplitude and the Hamiltonian density of the cw is shown for a range of magnitudes of the BEC fields 
with Mf = ±1. The (dimensionless) reference amplitude is ^ref = — [?i^(fci — fc_i)^/8m + qB^]/c2. 

coupled nonlinear Schrodinger equations, without any parametric terms that either mix the spin components or lock 
their frequencies (chemical potentials) together. If such solutions are viewed in a rotated reference frame that mixes 
the spin components, there is an oscillation due to the beating of one frequency and wavenumber against another. 
The oscillations with frequencies ([T5|) . in contrast, cannot be eliminated by a change in the variables. There is 
just one oscillatory frequency for a given near-cw solution. One could rotate such a ncar-CW solution with a small 
oscillation into another reference frame, and get oscillations as in Refs. [13, on top of the frequency (fT6|) of the 
small variations about the fixed point. The oscillations in Refs. (33 | can be affected by linear Zeeman splitting; the 
oscillation frequencies (|16p are independent of it. 

IV. SUMMARY AND CONCLUSIONS 

We obtained the most general continuous wave (plane wave) solutions for spinor BECs for spin F = 1, with linear 
and quadratic Zeeman splitting due to a magnetic field. We do not make the assumption that the wavenumbers of 
the different spin components are all the same. The physics only requires that cw solutions have wavenumbers and 
frequencies in the Alp = fields that are the average of the quantities in the Mp = ±1 fields. There are three 
distinct families of cw solutions. The first family are the solutions in which the Mp = spin component vanishes. 
The second and third families of cw solutions have non-zero densities of particles with spin Mp = 0. The second and 
third families of solutions are distinguished from each other by the phase of the Mp = BEC field with respect to the 
average phase of the Mp = ±1 BEC fields. Other things being equal, the different phases in the Alp = fields give 
the cw solutions different densities of spin Mp = particles. The third family does not allow all the wavenumbers to 
be identical; in this limit, the density of Mp = particles is asymptotically large. 

The first family of solutions, without spin Mp = particles, is governed by the coupled nonlinear Schrodinger 
equations, which have been studied in great detail. The solutions of the second and third type depend on whether 
the BEC is ferromagnetic or polar (antiferromagnetic). There are always cw solutions of the first type, and there is 
always at least one solution of the second or third type. 

The linear Zeeman splitting term, with a constant uniform magnetic field, can be eliminated by a gauge transfor- 
mation, so it does not change the dynamics. Quadratic Zeeman splitting cannot be eliminated, and has nontrivial 
effects. If the quadratic Zeeman coefficient is positive (which is the usual case for F = 1 particles), it's effects on the 
cw solutions are similar to the effects of a difference in the wavenumbers in the spin components. If the quadratic 
Zeeman coefficient is negative (which can be achieved by using microwaves and the Stark effect), a sufficiently large 
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magnetic field can made the cw solutions in the F = 1 spinor BEC switch from ferromagnetic to antiferromagnetic 
and vice versa. 

The cw solutions of the second and third families are stable to small perturbations within the cw ansatz. We 
calculated the frequency of these perturbations, which correspond to oscillatory coherent spin mixing. 
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